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Abstract 

In the search for tidal structure in Galactic satellite systems, we have conducted a photometric survey over a 
10 deg 2 area centred on the Fornax dSph galaxy. The survey was made in two colours, and the resulting colour- 
magnitude data were used as a mask to select candidate Fornax RGB stars, thereby increasing the contrast of 
Fornax stars to background sources in the outer regions. Previously, we reported the presence of a shell (age 
2 Gyr) located towards the centre of Fornax. In this contribution we reveal a second shell, significantly larger 
than the first, located 1.3° NW from the centre of Fornax, outside the nominal tidal radius. Moreover, the 
distribution of Fornax RGB stars reveals two lobes extending to the spatial limit of our survey, and aligned 
with the minor axis and with the two shells. These results support the hypothesis of a merger between Fornax 
and a gas-rich companion approximately 2 Gyr ago. 
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1 Introduction 

Observations of large-scale structure throughout the Uni- 
verse indicate galaxies are arranged in clusters and fila- 
ments, surrounding large voids which contain few luminous 
objects. The processes leading to this structure have been 
investigated by dark matter-dominated cosmological simu- 
lations (for examp le, iNavarro. Frenk. fc WhTtelll996l Il997l 
iMoore et al 1ll999h which predict that a complex system of 
hierarchical merging formed the galaxies of the present day. 
Evidence of these mergers can be seen a s remnant struc- 
ture in large spiral and elliptical galaxies. iMalin fc Carter! 
( 1980) detected phase-wrapped shells in the haloes of 
large gal axies throug h a process of photographic ampli- 
fication jMalinlll978li. and found that such structure is 
relatively common j Malin fc Carteilll983l) . Simulations by 
iHernauist fc Quinnl (ll98^^T989T) indicate that these shells 
originated through the merging of a smaller companion 
galaxy with its host object, with some of the companion's 
mass (~10%) distributed as shells moving on radial or- 
bits through the potential of the large galaxy. Hence shell 
structure is a common remnant of large galaxy mergers. 

If such structure can be seen in large galaxies, it is 
reasonable to expect a similar effect in smaller bodies such 
as the dwarf galaxies. CDM simulations predict that these 
objects have also formed through a process of hierarchical 
merging, and are therefore a compilation of several clumps 
of matter. Also, the process of galaxy-galaxy interactions 
are thought to be scale-free; hence the formation of shell 
structure at large galactic scales should be echoed at the 
scale of dwarf galaxies. However, by redshift z ~ 10 dwarf 



galaxies are expected to have formed through the merger 
of these small dark halos, and late infall in a typical dwarf 
halo is thought to be negligible (Power 2003, personal com- 
munication). Thus, a typical dwarf galaxy is not expected 
to display recently formed shell structure. 

Based on a deep photometric dataset covering 
the inner 1/3 deg 2 of the Fornax d Sph obtained by 
IStetson. Hesser. fc Smecker-Haneilll99^ . we previously re- 
ported the discovery of a shell-like feature located 17 ar- 
cmin (~1.8 core rad ii) from the centre of the system 
iColeman et alJ 120041 : hereafter C04). The analysis re- 
vealed that the clump is composed almost exclusively of a 
stellar population with an age of 2 Gyr. It is aligned with 
the major axis of Fornax, and is situated on the minor axis. 
Our inference was that the stellar association represents 
shell structure in a dwarf galaxy, the remnant of a merger 
between Fornax and a gas-rich companion approximately 2 
Gyr ago. However, such a conclusion was based on a single 
shell located towards the central regions of Fornax. 

In this contribution we report the discovery of a second 
shell, located outside the nominal tidal radius of Fornax, 
based on a wide- field photometric survey of the dSph. This 
shell is oriented parallel to the original structure described 
in C04, and is located on the opposite side of Fornax. In ad- 
dition, two large structures aligned with these shells extend 
approximately two tidal radii from the centre of Fornax. 
Based on these results, the system appears to have experi- 
enced a significant interaction with a gas-rich companion 
in the recent past. A full description of the Fornax survey 
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Figure 1. The colour-magnitude data from the inner 40' of For- 
nax, compiled from our photometry survey. The red lines out- 
line the selection region for candidate Fornax stars. The nominal 
photometric limit at V = 20 represents the limit of the entire 
datasct, while the deep photometric limit (V = 20.7) applies for 
six of the sixteen fields. 



and anal ysis techniques can be found in IColeman et alJ 
flin prep.|). 



Figure 2. Distribution of stars in the V20.7 dataset, which only 
includes six of the sixteen fields imaged. The full survey area is 
marked as the dotted line, where the remaining ten fields have 
been excised from this image. The points plotted are candidate 
RGB stars for the six fields with complete photometry down the 
the deep photometric limit displaye d in Fig . HI The red ellipses 
define the core and tidal radii from iMateol Jl99St) . The original 
shell described in C04 is located at Aa = 0.10°, A8 = -0.25° 



2 The Fornax Survey 

We have collected photometry in V and I bands for an 
area of 3.1 x 3.1 deg 2 centred on the Fornax dSph using 
the Wide-Field Imager (WFI) attached to the SSO 1 metre 
telescope. The region was imaged using a mosaic of 4 x 4 
fields (labelled Fl, . . . , F16), and the data are complete to 
a magnitude of V — 20 at the colour of the RGB. All fields 
contain an overlap region (width ~5') with each adjacent 
field, which we used to measure the internal accuracy of 
the photometry and astrometry. These analyses, as well as 
the coordinat es and a map of the survey region are pre- 
sented in lColeman et all din prep J) . The colour-magnitude 
diagram (CMD) from the inner 40' of Fornax is shown in 
Fig. where the upper dotted line illustrates the com- 
pleteness limit of the survey. By selecting stars within the 
RGB selection range, we are able to increase the 'signal' 
of Fornax stars compared to the background field popula- 
tion (the 'noise'). We selected candidate RGB stars down 
to the completeness limit, and label this the 'V20' dataset. 
However, six of the sixteen fields we imaged were com- 
plete to a deeper limit of V = 20.7, represented as the 
lower dotted line in Fig. Hence we selected stars down 
to this deeper limit in these six fields to further increase 
the signal-to- noise of Fornax stars, and label this the 'V20.7' 
dataset. A full discussion of both datasets is presented in 
IColeman et alJ ilin prep.|) . 

In Fig. |21 we show the distribution of the stars in 
the V20.7 dataset. The outer red ellipse de fines the nom - 
inal tidal radius of this system, taken from Mateo (1998). 



Perhaps the most apparent feature in Fig. [21 is a region 
of increased stellar density located approximately 1.3° 
north-west from the centre of Fornax. This apparent over- 
density is situated in the field F4 outside the nominal 
tidal radius. An inspection of the CCD images for F4 
indicates all the stars within this feature are real stellar 
sources, and the star finding algorithm does not appear 
to have missed any stars in the remainder of the field. 
Also, given the high Galactic latitude of Fornax, we do 
not expect this feature to be caused by dust extinction. 
Corre spondingly, a visual inspection of th e IRAS 100/im 
map (Schlcael. Finkbeiner. fc DavisHl99al does not show 
any significant dust structure in the region of this feature. 

Thus, the feature merits further analysis. We have 
measured the total surface density in the vicinity of the 
feature to be 1.7 ± 0.1 times that of the remainder of 
F4, where the region we have selected is outlined in Fig. 
|3] That is, at least one third of the stars in the over- 
dense region are part of the feature, while the remainder 
form the background popul ation. We remov ed the back- 
ground population IColeman et alJ ilin prep.h and, assum- 
ing the feature lies at the same d istance modulus as Fornax 
[(m — M) = 20.70: lMateolll99S|| . then it has an integrated 
absolute magnitude of Mv = — 6.7±0.2, comparable to the 
five globular clusters associated with Fornax. A visual in- 
spection of_Hie_£olovrr^agjrh^ude data for all stars in this 
feature iColeman et alilin prepJ) reveals that the stellar 
population defines an RGB qualitatively matching that of 
Fornax. 
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Figure 3. Distribution of stars in the V20.7 dataset for the NW 
quadrant only. The red box outlines the shell selection region. 



3 Analysis 

A full description of the an alysis techniques are presented 
in IColeman et aD llin prep J) . Here we describe one of the 
tools used to statistically investigate the extra-tidal stellar 
structure. We constructed a Monte Carlo simulation where 
each iteration consisted of placing a circle (radius 12 ar- 
cmin) at a random position outside the nominal tidal ra- 
dius and counting the number of stars therein. By complet- 
ing 10 6 iterations, we were able to construct a function de- 
scribing the probability of measuring a given stellar density 
beyond the Fornax tidal limit. In order to completely sam- 
ple the survey area, this simulation was conducted over the 
V20 dataset. The star counts within circles falling partially 
outside the survey region were scaled by the appropriate 
amount to correct for the non-surveyed area. The resulting 
probability function is displayed in Fig. [I] A random distri- 
bution of stars would be expected to display a normal prob- 
ability distribution in density. However, immediately visi- 
ble in Fig. g| is a triple-peaked structure, thus implying the 
extra-tidal region consists of three populations with distin- 
guishable densities. The dotted line corresponds to the pre- 
dicted density of the field sta r populat ion, determined by 
iRatnatunga fc Bahcall <| l985|) u sing the[B ahcall fc Soneiral 
( 1980) Galactic model. IRatnatunga fc Bahcalll jl985l) esti- 
mate the star count uncertainty at 25%, hence we conclude 
the low-density peak corresponds to the true field star den- 
sity. We measured the standard deviation from the mean 
density of the field population, and this value (erf) allows us 
to measure the statistical significance of the higher density 
structures. 

The high-density peak (centred at a density of 1150 
stars arcmin" 2 ) corresponds to the feature described in 
Sec. 2, situated at Aa = -1.0°, AS = 0.8° in Fig. H This 
appears to be the only high-density structure beyond the 
nominal tidal radius of Fornax. With a mean stellar density 
~20(jf higher than that of the field population, the feature 
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Figure 4. Monte Carlo simulation of the density distribution of 
Fornax RGB-selected stars, showing the probability of measur- 
ing a given density beyond the nominal tidal ra dius. The d ashed 
line c orresponds to the (normalised) IRatnatunga fc Ba hcall 
1 1985) star counts. The red line indicates the sum of two best-fit 
Gaussians to the low and medium density background popula- 
tions. 

is clearly separated from the background. It is located on 
the opposite side of Fornax to the previous shell (C04), and 
both are situated on the minor axis and aligned parallel to 
the major axis. Therefore this structure appears to be a 
second Fornax shell. 

Finally, we examine the medium-density population. 
The contour plot in Fig. [^1 was constructed by convolving 
each star in the V20 dataset with a Gaussian of approxi- 
mate width 10'. The contour smoothing length is 3 arcmin. 
The lower significance limit is 4.5o"f , hence the first contour 
traces all structure with a density 4.5o"f above that of the 
field star population. This contour outlines two 'lobes' ex- 
tending beyond the nominal tidal radius of Fornax, located 
on the minor axis. Assuming they are at the same distance 
as Fornax, with a similar stellar population, they measure 
~8 kpc in length with an integrated absolute visual mag- 
nitude of Mv ~ — 9. 

In summary, beyond the nominal tidal radius of the 
dSph, the NE/SW quadrants of the survey contain only 
the field star population, with no discernible excess of For- 
nax RGB stars. However, in the NW/SE quadrants, we 
detected two large structures extending at least 1.0 deg be- 
yond the tidal limit, and these appear symmetrical along 
the minor axis of Fornax. Moreover, we have located a sec- 
ond shell approximately 1.3° NW of the centre of Fornax. 



4 Discussion 

In C04 we proposed the existence of shell structure in For- 
nax based on a single feature located approximately two 
core radii from the centre of the system. Such structure 
would arise from the interaction between Fornax and a 



© Astronomical Society of Australia 0000 



4 M. G. Coleman & G. S. Da Costa 



T 




J i i i i I 

1 

Ac C) 



Figure 5. Distribution of Fornax RGB stars in the V20 datasct 
where each star has been convolved with a Gaussian of approx- 
imate width 10'. The contours are logarithmically spaced, and 
are fitted with a smoothing length of 3.0'. The first contour rep- 
resents a stellar density 4.5<Tf above the field star population. 



companion object. Another possibility is that the two large 
lobes of material extending to either side of the dSph rep- 
resent tida l tails. They qualitatively match simulat ions of 
tidal tails jHelmi fc Whitell200lt[Maver et al]l200ll) which 
are caused by interaction with the Milky Way. These mod- 
els also indicate that material leaves the satellite in bursts, 
which may produce the observed shell-like structure. 

However, clear tidal tails have only been observed in 
a single dSph satellite of the Galaxy; the Sagittarius dSph. 
The Sagittarius orbit brings it within approximately 15 
kpc from th e Galactic centre with an orbital e ccentricity 
of e ~ 0.75 llLaw. Johnston, fc MaiewskilfeoO^ . and thus 
Sgr experiences substantial tidal forces. In contrast, For- 
nax is located at a distance of 140 kpc (lMated ll998') and 
proper motion measurements indicate it is on a polar or- 
bit with a small eccentricity (e = 0.2 7) and is currently 
at perigalacticon iDinescu et al.ll2004l) . Moreover, Fornax 
contains a mass larger than all the other Galactic dSphs 
combined (excluding Sagittarius), which further precludes 
the possibility of stars escaping from this relatively deep 
potential well. Thus, given the relatively large mass, dis- 
tance and orbital parameters of Fornax, we do not expect 
to observe such clear evidence of Milky Way interaction. 

Further, the inner shell described in C04 demon- 
strated a clear signature; it was dominated by stellar 
population with an age of 2 Gyr, significantly differ- 
ent to the mixed stellar population ob served in Fornax 
iStetson. Hesser. fc Smecker-Hanelll99ct) . If the observed 
substructure is due to the gravitational influence of the 
Galaxy, then it would contain the same stellar population 
as the remainder of the system. 

The high density structure located 1.3° NW of Fornax 
is clearly shell-like in appearance. In the case of shell struc- 



ture, the lobes would represent those portions of shell for 
which we do not see a limb brightening effect; in essence, 
they lie in front of, or behind, the Fornax system. These 
appear remarkably similar t o the simulations of interact- 
ing dark halos presented bv iKnebe et alJ ll2004h . To fur- 
ther investigate the shell possibility, we will obtain deep 
photometry of the shell region. If the stellar population is 
similar to the Fornax system, this would indicate the sub- 
structures are tidal tails. However, if the stars in the shell 
display a marked age difference from those in Fornax, this 
would support the possibility of shell structure in a dSph. 
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